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The recent introduction of the concept of scattering fidelity, causes us to revisit the experiment 
by Lobkis and Weaver [Phys. Rev. Lett. 90, 254302 (2003)]. There, the "distortion" of the coda 
of an acoustic signal is measured under temperature changes. This quantity is in fact the negative 
logarithm of scattering fidelity. We re-analyse their experimental data for two samples, and we 
find good agreement with random matrix predictions for the standard fidelity. Usually, one may 
expect such an agreement for chaotic systems only. While the first sample, may indeed be assumed 
chaotic, for the second sample, a perfect cuboid, such an agreement is more surprising. For the first 
sample, the random matrix analysis yields a perturbation strength compatible with semiclassical 
predictions. For the cuboid the measured perturbation strength is much larger than expected, but 
with the fitted values for this strength, the experimental data are well reproduced. 



Lobkis and Weaver (henceforth LW) experiment 
have measured the sensitivity of elastic coda waves to 
temperature changes. After correcting for a trivial term 
due to change of volume and wavespeed they quantify 
the effect as "distortion" , and study its decay as a func- 
tion of time. We shall show that this quantity is the 
negative logarithm of the "scattering fidelity" introduced 
and measured in 0, 13 ■ For sufficiently chaotic dynamics 
in systems weakly coupled to decay channels, scattering 
fidelity approaches the standard fidelity amplitude. Fi- 
delity, which is the absolute value squared of the fidelity 
amplitude, has received a great deal of attention in recent 
years, as it is used as a benchmark in quantum informa- 
tion processes 0, and in the context of quantum chaos 
(for a partial overview and references see |5|)- 

For the above mentioned reasons, we expect the ran- 
dom matrix analysis of fidelity decay , to apply to the 
experiments of LW, as it does in the case of comparable 
experiments with chaotic microwave cavities |2, Q • Ran- 
dom matrix theory (RMT) makes a definite statement on 
the form of the fidelity amplitude as a function of time. 
It yields a unified description of the "perturbative" and 
the "Fermi golden rule" 0, regime. This analysis is 
more appropriate than the treatment in Q , which seems 
to be on a footing similar to the Fermi golden rule re- 
sults obtained in . RMT still contains the perturbation 
strength as a free parameter, which must be determined 
independently. For quantum chaotic systems, different 
semiclassical methods have been used 8, 9]. It is still 
an open problem to extend those results to the pertur- 
bative regime Q- 

We analyzed the data for two samples, measured by 
LW, both being three dimensional objects of similar size 
in all three spatial directions. The first object is the so 



called "medium block" (after extra cut) which is sup- 
posed to have dominantly chaotic dynamics, with no 
symmetries left. The second object is a cuboid, called 
"rectangle", where the dynamics is not chaotic, but due 
to mode conversion it is also not integrable and actu- 
ally known to display random matrix behavior |lOl |. In 
both cases, wc obtained very good agreement with our 
RMT result, as far as the form of the scattering fidelity 
functions is concerned. For the first sample, we also ob- 
tain a reasonable agreement for the perturbation strength 
(comparing the fitted values from an RMT analysis with 
the theoretical result of LW). For the second sample, we 
obtain values which are systematically too large, by a fac- 
tor of two, approximately. Over all, we find that our ran- 
dom matrix approach provides a very good framework, 
for the discussion of distortion measurements of non- 
integrable elastodynamic systems, thus providing addi- 
tional information for the elastic problem. Conversely, 
experimental results with adequate control of the per- 
turbation strength are rare in fidelity analysis and thus 
the possibility to use acoustic scattering results is highly 
welcome. 

Scattering fidelity The scattering fidelity Q is defined 

as 

fab{t) = {S,,{tr SUt))/^{\Sa,{t)?){\SU{t)?) . (1) 

Here, {Sab{t)* S'^,^{t)) is the Fourier transform of the 
cross-correlation function of a scattering matrix element 
for the unperturbed and perturbed system, respectively. 
Similarly, (|S'afc(t)P) and {\S'^^(t)\'^) are the Fourier trans- 
forms of the corresponding autocorrelation functions, 
used for proper normalization. It can be shown that 
appropriate averaging (denoted by (. . .)) yields the stan- 
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dard fidelity amplitude for chaotic systems weakly cou- 
pled to the scattering channels 0- In that RMT 
model can be used to compute the fidelity amplitude f{t) 
on the basis of linear response theory @. With x = t/tn 
(tn is the Heisenberg time) and ^2(2;), the two point form 
factor for the Gaussian orthogonal ensemble, one obtains: 



dx' dx" b2{x") 

Jq 



(2) 



Here, the "perturbative" and the "Fermi golden rule" 
regime arise as two particular limits: For small Aq the 
decay of f{t) is dominantly Gaussian, while for large Aq 
it is dominantly exponential 0. In a microwave exper- 
iment excellent agreement between fab(t) and f{t) has 
been found HQ. The validity of Eq. © for the fidelity 
of closed chaotic systems has been demonstrated with 
numerical calculations for the kicked rotor • 

In the authors measure the acoustic response to 
a short piezoelectric pulse as a function of time. They 
consider the normalized cross correlation between two 
such signals obtained at different temperatures Ti and 
T2: 



X{e) 



J dt SrAt) StM^ + e)) 
^JdtS^^^it) JdtSl^{t{l + e)) 



(3) 



This equation displays a very similar structure as Eq. 
The time averaging over a small window, performed in 
Eq. 0, corresponds to a smoothing of the correlation 
functions in Eq. The selection of e, such that the 
correlation function X{e) ~ ^max becomes maximal, 
eliminates the trivial effects due to dilation and change 
of wavespeed, caused by the temperature change. This 
is quite similar to the spectral unfolding performed in 
Ref. to eliminate the volume effect of a moving wall. 
The distortion is defined as D{t) = — ln(Xijiax)j where 
the time dependence is given by the "age" of the signal, 
i.e. the center of the small time- interval over which the 
correlation function X(e) was evaluated. If formulated as 
a scattering process, we find D{t) = — ln[/aa(t)], where 
the scattering channel a is defined by the transducer, that 
transmits excitation to and from the sample. Thus, for 
sufficiently chaotic samples in the elastodynamic scatter- 
ing experiments, we expect that the scattering fidelity is 
equal to the fidelity amplitude, and that the latter is well 
described by the RMT result, Eq. 0. 

In Q the perturbation strength is estimated in 
terms of the distortion coefficient C « 3.26 x 
10~^°/(K^ cmmsec""'^). This value is based on a ray pic- 
ture of the resonating acoustic waves, the assumption of 
random reflection angles along the ray paths, and an es- 
timate for the mode conversion rates between dilational 
and shear waves. We use that value to obtain an estimate 
for the dimcnsionlcss perturbation strength Aq in Eq.(|51l: 
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FIG. 1: (Color online) The distortion for the medium block as 
plotted in Fig. 4 of Ref. The thin jagged lines correspond 
to measurements in the frequency ranges 300 kHz, 600 kHz, 
700kHz, and 800kHz (from bottom to top). The thick lines 
show — ln/(t) according to Eq. J^J, using the values of Ao, 
as given in Tab. 1 (full RMT fit). For the thin straight lines, 
only the linear term in Eq. Q has been taken into account. 



where A denotes the temperature difference Ti — T2 in 
Kelvin, while V and S denote the volume and the sur- 
face of the sample. From Refs. [mil, we obtain the 
following expression for the Heisenberg time tn- 
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Ao = / V^ctnif) c = C A^ V/S 



(4) 



Here, the longitudinal wave velocity is Cp = 637 cm/msec 
and the transverse (shear) wave velocity is Cs = cp/q = 
316 cm/msec. The frequency (range) is denoted by /. 

Medium block That specimen has volume V — 
906 cm'^ and S = 636 cm^, and the temperature differ- 
ence was A = 4K. The distortion, measured as a function 
of time is shown in Fig. 1. With the Heisenberg times 
given in Eq. |SJl, we may fit the perturbation strength Aq 
and compare to the estimate in Eq. (0J. For experimen- 
tal reasons, the data for D{t) is reliable for t > 20 msec, 
only. The fits, mentioned above, have been restricted 
correspondingly. 

The agreement between the measured distortions or 
scattering fidelities and random matrix theory is within 
the statistical error of the data. A deviation from lin- 
ear behavior is hardly noticeable, except for the 300 kHz 
data, which is the only case where tn « 78 msec ^3 lies 
within the time range of the figure. However, there is still 
a considerable effect on the fit values for Aq. This can be 
seen from the dashed curves, which have been obtained 
by taking only the linear part of Eq. (2) into account. 

Tab. 1 gives the values of Aq as obtained from linear 
fits to the data (such a procedure has been applied in Q), 
and from fits with the full random matrix result, Eq. 
In the last column, the respective Heisenberg times are 
given. In Fig. 2, the fitted perturbation strengths are 
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//kHz 


Ao (lin. fit) (full fit) iff /msec 


300 


0.471 


0.298 


78.4 


600 


1.584 


1.381 


294.6 


700 


2.136 


1.929 


397.3 


800 


2.709 


2.505 


515.3 



TABLE I; Medium block: Table of linear and full RMT fit 
values for the dimensionless perturbation strength Ao together 
with the respective Heisenberg times for the four frequency 
ranges, analysed. 
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FIG. 2: Medium block: The perturbation strength as a func- 
tion of the frequency range. The circles show the values for 
Ao obtained from a fit with Eq. (|5J , while the diamonds show 
the corresponding values when only the linear term of Eq. Q 
is taken into account. The solid line gives the perturbation 
strength as obtained from Eq. 

plotted versus frequency. The solid line shows the esti- 
mated behavior of the perturbation strength, computed 
from Eq. As long as the volume term dominates 

in the expression Q for tn, the perturbation strength 
increases quadratically with the frequency. The differ- 
ence between the linear fits (diamonds) and the full RMT 
fits (circles) is clearly noticeable. The full RMT analysis 
moves the fit values for Xq quite close to the prediction 
of LW. We have no explanation for the remaining de- 
viations. It is not clear, whether these are statistically 
acceptable, whether chaoticity is not perfect, or wether 
there is some other reason. We may recall that in Ref . Q , 
the experimental perturbation strength did also not agree 
with the theoretical estimate, as it was too small. How- 
ever, there, the perturbation strength could be measured 
independently, from the level dynamics. This demon- 
strated that for the cavities under study, the semiclas- 
sical approximation was not yet justified. It would be 
desirable, if the level dynamics could be measured in the 
context of the present experiment. 

Rectangle In Fig. 3 and Fig. 4, we analyse data for 
the rectangle, a perfect cuboid. Clearly, a scalar wave 
equation would lead to integrable ray dynamics, where 
our random matrix model must fail. However, in the 



//kHz 


Ao (lin. fit) (full fit) tn/msec 


100 


0.241 


0.061 


10.36 


200 


0.607 


0.279 


35.75 


300 


0.903 


0.567 


76.19 


400 


1.389 


1.033 


131.68 


500 


1.905 


1.561 


202.21 


600 


2.564 


2.227 


287.78 


700 


3.282 


2.957 


388.39 


800 


3.969 


3.662 


504.04 



TABLE II: Rectangle: Table of linear and full RMT fit values 
for the dimensionless perturbation strength Ao together with 
the respective Heisenberg times. 
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FIG. 3: (Color online) The distortion as a function of time, 
for the rectangle. The thin jagged lines correspond to mea- 
surements in the frequency ranges from 100 kHz to 800 kHz, 
in steps of 100 kHz (from bottom to top). The thick smooth 
lines show the best fits with — ln/(t), according to Eq. (2), 
with Ao given in Tab. 2 (full RMT fit). 

present case, the wave field has two components (dilata- 
tional and a shear waves) , which are coupled due to mode 
conversion. The corresponding classical dynamics are 
marginally stable, but may still be ergodic. Recently, 
the elastodynamic system has been studied thoroughly, 
with a focus on spectral statistics Apart from cer- 
tain symmetries, the statistical measures show accurate 
RMT behavior. Nevertheless, from a theoretical point 
of view, it remains an open question, whether scattering 
fidelity still agrees with standard fidelity, and whether 
the behavior of the fidelity amplitude can be described 
by random matrix theory. The following analysis with 
Eq. 121) will shed some light on these questions. 

For the rectangle, we have data for eight frequency 
windows at / = 100 kHz, 200kHz, 800kHz. The 
corresponding Heisenberg times range from 10 msec to 
500msec (see Tab. 2). In Fig. 3, one can clearly see a 
transition from linear to quadratic decay, characteristic 
for the RMT expression, Eq. (j2Jl. Here, it is really sur- 
prising that the RMT expression describes the data so 
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FIG. 4: Rectangle: The perturbation strength as a function 
of the frequency range. The circles show the values for Ao as 
obtained from Eq. and from the same expression, taking 
into account only the hnear term (diamonds). The solid line 
gives the perturbation strength as given by Eq. @. The 
dashed line shows the same expression, but with C = 8.90 x 
10~^''/(K^ cmmsec"^) obtained from a fit to the data points 
(circles) . 

well. The perturbation strengths obtained from fits to 
the RMT expression are plotted in Fig. 4 (circles). Ex- 
cept for a constant factor (the distortion coefficient), the 
result follows the theoretical expectation (Eq. solid 
line). A fit for the distortion coefficient on the basis of 
Eq. I0J yields the dashed line. The values for Aq ob- 
tained from a linear fit (diamonds) clearly, do not follow 
the theory. 



In this work, we identified a previously published elas- 
tomechanic scattering experiment as an experiment that 
measures scattering fidelity in a setting usualy called echo 
dynamics. As in the case of electromagnetic billiards, it 
turns out that a simple random matrix model describes 
the decay of the scattering fidelity very well, despite of 
the fact that these are much more complicted systems. A 
number of questions could not be answered. In the case 
of the medium block, we had not enough data, in the 
low frequency region, where the crossover between Fermi 
golden rule and perturbative regime occurs. Also, we 
could not consider any possible symmetries of the system, 
which might be one source for the remaining differences 
between estimated and measured perturbation stren gths . 
For the rectangle, our results complement those of |10l| . 
They show that at least in one important aspect the wave 
functions behave like those of a chaotic system. 

The high quality factors of elastomechanic experi- 
ments, as well as the possibility to measure cxplicitely 
in the time domain, make these experiments particu- 
larly atractive. Among possible experiments, we believe 
that it will be worthwhile to analyze strong perturbation 
data, e.g. larger temperature differences. Such experi- 
ments are more difficult, due to precision problems, but 
may serve to explore the limits of RMT, or detect other 
regimes of fidelity decay. This is particularly interest- 
ing, because the exact solution of the RMT model is now 
available [T5| . 
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